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We describe a family of icosahedral structures for lead clusters. In general, structures in this 
family contain a Mackay icosahedral core with a reconstructed two-shell outer-layer. This family 
includes the anti-Mackay icosahedra, which have have a Mackay icosahedral core but with most 
of the surface atoms in hexagonal close-packed positions. Using a many-body glue potential for 
lead, we identify two icosahedral structures in this family which have the lowest energies of any 
known structure in the size range from 900 to 15000 lead atoms. We show that these structures are 
stabilized by a feature of the many-body glue part of the interatomic potential. 



I. INTRODUCTION 

The structure of an atomic cluster often differs from 
that of the corresponding bulk material]. In such a clus- 
ter, the number of surface atoms is comparable to the 
number of interior atoms, and consequently, the surface 
energy plays an important role in determining the overall 
structure. For example, regular noncrystalline structures 
with fivefold axes of symmetry, such as icosahedra and 
decahedra, are known to occur in gold .and a variety of 
other face-centered cubic (fee) metalsErBfl Such struc- 
tures are comprised of deformed fee tetrahedral units 
where adjacent tetrahedral faces meet at a twin plane. 
The energy cost of twinning at the inner tetrahedral 
faces, and the strain energy in the deformed tetrahedra, 
is overcome by the-resulting energetically favorable close- 
packed outer facesotl 

The structure of a cluster, while not only of fundamen- 
tal interest, is also a key determinant of many of its prop- 
erties. However the delicate balance between surface and 
internal energies often produces a complex dependence 
of structure upon cluster sizeall. Eventually, as the size 
of a cluster increases, the bulk structure must win out, 
but at sizes below this, clusters can assume a variety of 
regular non-crystalline structures. For systems that can 
be adequately described by pair potentials, there is a rel- 
atively good understanding of how structure depends on 
the form of the potentiaEiJ. However, for many systems 
of interest, such as metals, the interatomic interactions 
are more complex. 

Metals exhibit a strong many-body character to their 
bonding, and because of this, the competition between 
fee, icosahedral and decahedral structures is less well un- 
derstood. In addition, many-body effects can potentially 
lead to the emergence of new structural formstJuU. In 
order to study the effect of this bonding on structure, it 
is often necessary to resort to empirical potentials, as ab 
initio electronic structure methods are prohibitively ex- 
pensive for all but the smallest cluster sizes. For exam- 
ple, lead clustess-have been studied using a many-body 
glue potentiall3't3. The hrst comprehensive computa- 



tional study of lead cluster structure, by Lim, Ong and 
Ercolessili-J, utilized this potential to compare the ener- 
getics of closed-shell fee cuboctahedral and Mackay icosa- 
hedral structures for cluster sizes from 55 to 3871 atoms. 
By directly comparing the binding energies at the same 
numbers of atoms for each structure, they demonstrated 
that the fee cuboctahedra were favored over Mackay 
icosahedra over this size range. 



However, electron diffraction of 3-7 nm lead clusters 
has produced diffraction patterns that|-eannot be ade- 
quately fitted by fee structural modelala. In addition, 
recent simulations of the melting and freezing of clus- 
ters, usipg the same glue potential as Lim, Ong and 
Ercolessilij, have unexpectedly revealed that fee struc- 
tures were not the lowest in energy in this size rangetS. In 
these simulations, clusters were prepared by the quench- 
ing of liquid lead droplets. This procedure was found 
to produce icosahedra overwhelmingly, and these resolid- 
ified icosahedra were, found to be energetically favored 
over fee structuresta. The resolidified icosahedra re- 
sembled anti-Mackay icosahedraliJ which have a Mackay 
icosahedral core but with most of the outer layer in 
hexagonal close-packed (hep) surface sites. The improved 
stability of these icosahedra was evidently due to this sur- 
face reconstruction. 



The purpose of this work is to examine the surface re- 
construction of these resolidified icosahedra in more de- 
tail. In particular, we will show that the surface recon- 
struction is similar, but not identical, to that of the anti- 
Mackay icosahedra. These surfaces features lead us to 
identify a new family of icosahedral structures, of which 
the anti-Mackay icosahedron is a member. We will de- 
scribe these structures in detail, and show how they can 
lead to lower energy structures for the lead glue potential. 
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II. COMPUTATIONAL TECHNIQUES 
A. Molecular dynamics 

Molecular dynamics was performed using a local ver- 
sion of the classical molecular dynamics code ALCMD, 
originally developed by Ames Laboratory. Finite tem- 
perature simulations were performed in the microcanon- 
ical ensemble (i.e. constant energy). The time step was 
chosen as 3.75 fs throughout. Melting and freezing sim- 
ulations were caccied out using the procedure detailed in 
Hendy and HalM 

The inter-atomic potential used is due to Lim, Ong 
and ErcolessitJ. This is a many-body glue-type potential, 
given by 

E = -Epair + E g i uc , 

= Z>(r«)+£tf(r*) , (1) 

i<j i 

where cj> is a short-range pair potential and U(n) is a 
many-body glue term which reflects the effects of non- 
local metallic bonding. The quantity rii is a "generalized 
coordination number" for atom i defined as: 

ni = ^2p{rij) (2) 

3 

where p is some short-ranged "atomic density" function. 
In practice, the function p(rtj) has a cut-off r cut , beyond 
which p(rij) = 0, and for the potential here, a value 
of r cut = 5. 503 A is used. This cut-off typically lies be- 
tween the second and third neighbor shells. The three 
functions (j), p and U have been obtained by fitting to 
a number of known properties of lead including cohe- 
sive energy, surface energy, elastic constants, phonon fre- 
quencies, thermal expansion and melting temperatureEl 
This po tantial has been used previously to model lead 
clustersESIlalla, temperature-dependent surface recon- 
structions of low-index lead surfaceslij and pre-melting 
of low-index lead surfacesel Recently, putative global 
minimum energy cluster structures were determined for 
this potential, for cluster sizes of up to 160 atomsei 

B. Common neighbor analysis 

Common neighbor analysisEiLpCNA) has been used 
here to analyze cluster structuresoE3. CNA is a decompo- 
sition of the radial distribution function (RDF) according 
to the local environment of each pair. We consider that 
the first peak in the RDF represents "bonded" neighbors. 
As such, if r c is the first minimum in the RDF, we classify 
any pair separated by r < r c as a bonded pair. With this 
identification, any pair contributing to the RDF can be 
classified by a set of three indices, ijk, which provide in- 
formation on the local environment of the pair. The first 
index, i, is the number of bonded neighbors common to 




FIG. 1: The anti-Mackay surface termination (left) has sur- 
face atoms in hep positions. The surface terminations con- 
sidered here do not include the vertex atoms (indicated by 
the dotted lines). The Mackay surface termination (right) is 
shown for comparison. 

both atoms. The second index, j, is the number of bonds 
between this set of common neighbors. The third index, 
k, is the number of bonds in the longest continuous chain 
formed by the j bonds between common neighbors. 

CNA is useful here because it allows one to distin- 
guish between local atomic arrangements, including fee 
and icosahedral environments, using the type and num- 
ber of ijk indices of each atom. For each atom, we can 
define riijk to be the number of bonds of this atom with 
CNA indices ijk. We can then classify the local envi- 
ronment of each atom using these n^k values. In Ta- 
ble | we have listed the classifications of CNA signatures 
used here to label the local environment of an atom (this 
classification is similar but not identical to that used by 
Cleveland, Luedtke and Landmann3). We note that these 
signatures are based only on the CNA decomposition of 
the first peak in the RDF. 

III. ANTI-MACKAY ICOSAHEDRA 

We will refer to an icosahedron with an anti-Mackay 
surface termination as an anti-Mackay icosahedron. With 
such a surface termination, atoms in the exterior shell lie 
in hep positions relative to a core Mackay icosahedron, 
as illustrated in figure [j] Note that for the terminations 
we will consider here, we neglect the icosahedral vertex 
atoms on the surface since these are not present in the 
resolidified icosahedra, and tend to increase the overall 
energy of the cluster. 

An n-shell anti-Mackay icosahedron contains a core 
(n — l)-shell Mackay icosahedron with l/2(n— 3)(n — 4) + 
3(n — 2) surface atoms per face (neglecting the 12 ver- 
tices). The total number of atoms in an n-shell Mackay 
icosahedron is 

Ico(n) = ^n 3 + 5n 2 + ^—n + 1, (3) 

giving the sequence 55, 147, 309, 561, 923, 1415, 2057, 
2869, 3871 ... Hence, the total number of atoms in an 
n-shell anti-Mackay icosahedron is 

Anti(n) = ^-n 3 + 5n 2 - ^-n - 1 (4) 
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TABLE I: Description of CNA signatures used. 



label 


description 


position 


classification of pairs 


A 


fee internal atom 


internal 


"■421 > 4 


B 


fee {111} face atom 


surface 


n3n > 3 and 71322 = 


C 


fee {100} face atom 


surface 


71211 > 3 


D 


fee {111}/{100} edge atom 


surface 


ri2ii = 2 and 71311 = 2 


E 


internal atom at a {111} fee stacking fault 


internal 


71422 > 5 


F 


internal icosahedral atom (spine or central atom) 


internal 


71555 > 2 


G 


surface icosahedral apex atom 


surface 


71555 > 1 


H 


surface icosahedral {111}/{111} edge atom 


surface 


71311 > 3, 71322 > 1 and 


1 


anti-Mackay surface {111}/{111} edge atom 


surface 


71311 = 2, 71200 = 2 and 


? 


unclassified signature (possibly disordered) 


internal /surface 





This gives a sequence of closed-shell anti-Mackay icosa- 
hedra with numbers of atoms as follows: 115, 267, 509, 
861, 1343, 1975, 2777, 3769, 4971, 6403, 8085 . . . 

Thus, the n-shell anti-Mackay icosahedron contains 
lOn+2 fewer surface atoms than the n-shell Mackay icosa- 
hedron, so the packing of the surface atoms will be less 
dense. However, these missing atoms come from edges 
and vertices with low coordination m . Table [il] compares 
the energy of an 8-shell Mackay icosahedron and an 8- 
shell anti-Mackay icosahedron, broken down by CNA la- 
bel. Here we can see that while the surface binding en- 
ergy per atom is worse for the anti-Mackay icosahedron, 
the binding energy per atom for the cluster as a whole 
is better than the Mackay icosahedron. Thus for the po- 
tential (|l|), this removal of edges and vertices from the 
Mackay icosahedron (which can be seen to have partic- 
ularly poor energetics in Table improves the total 
energy per atom. 

However, despite this improvement in the binding en- 
ergy of the anti-Mackay icosahedra, they are still not en- 
ergetically competitive with the cuboctahedra sequences 
or the resolidified icosahedra, as shown in figure g (note 
that here we refer to cuboctahedra with triangular (111)- 
faces simply as cuboctahedra, and to cuboctahedra with 
hexagonal (lll)-faces as truncated octahedra). Thus, 
while the resolidified icosahedra and the anti-Mackay 
icosahedra share a similar surface reconstruction, the 
resolidified icosahedra have other features which account 
for their more favorable energetics. We will discuss these 
features in the next section. 



IV. RESOLIDIFIED ICOSAHEDRA 

Hendy and HallB conducted a series of resolidification 
trials, where lead clusters were melted and then resolidi- 
fied at constant energy. Figure || shows the distribution 
of binding energies of 2057-atom clusters that emerged 
from a typical sequence of 25 resolidification trials (more 
details can be found in Hendy and Hallt3). These trials 
typically produce icosahedron-like structures, similar to 



0.020 r 



Mackay 
Icosahedra 




FIG. 2: Energies of clusters versus size: anti-Mackay icosa- 
hedra (solid square), cuboctahedra (solid diamond), Mackay 
icosahedra (solid triangle) and truncated octahedra (solid 
sphere). The energies are given relative to a fit to the en- 
ergies of the cuboctahedra sequence: i? cu boct — — 2.0293A + 
1.8216A 2/3 +0.7134iV 1/3 . 



that shown in figure |], which have higher binding energies 
than comparably-sized fee clusters such as cuboctahedra 
or truncated octahedra. In fact, the cluster shown in 
figure |] had the highesi-hinding energy produced in the 
trial. Hendy and HalO noted that these icosahedron- 
like clusters had similar surface reconstructions to the 
anti-Mackay clusters. 

In table O, we compare the structures of an anti- 
Mackay icosahedron and a resolidified icosahedron of a 
similar size. A crucial difference is the denser surface 
packing of the resolidified icosahedron. The number of 
atoms on the surface is comparable to that of the Mackay 
icosahedron (see table |J) , without the inclusion of ener- 
getically unfavorable edge atoms. Instead, there are a 
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TABLE II: Comparison of average energies of a closed-shell Mackay icosahedron, a closed-shell anti-Mackay icosahedron, a 
resolidified icosahedron and a m = 7 surface-reconstructed icosahedron. 



Cluster 







Mackay 


Anti-Mackay 




Resolidified 




m — ( 






Icosahedra 




Icosahedra 




Icosahedra 




Icosahedra 


Signature 


IN 


\E) 


N 


(E) 


IN 


\E) 


tvt 
IN 


\E) 






(eV/ atom) 




(eV/atom) 




(eV/atom) 




(eV/atom) 


A (tec) 


700 


-2.0159 


400 


-2.0264 


468 


O f\f\ A O 

-2.0048 


467 


-2.0062 


"D h 1 1 \ 

B (111) 


420 


-1.6254 


200 


-1.5836 


173 


-1.6963 


185 


-1.6780 


C (100) 






1 


-1.4007 


75 


-1.6425 


82 


-1.6254 


D (fee edge) 






- 


- 


2 


-1.2624 






li (ncp) 


fiOA 
630 


o m no 
-2.U103 


750 


-2.0071 


a c a 
654 


-z.U0o2 


046 


-2.U083 


F (ico) 


o c 
00 


-1.9916 


265 


-1.9794 


1 OQ 

loo 


-2.0016 


195 


i noon 
-1.9990 


Vjr IbLlIl ICO 1 


1 9 








1 

LO 


1 78^9 


i 


1 7709 


H (ico edge) 


210 


-1.4633 






1 


-1.9126 






I (anti edge) 






298 


-1.5531 


230 


-1.6301 


276 


-1.6212 


? surface 






61 


-1.3585 


137 


-1.5031 


83 


-1.6652 


? interior 










119 


-1.9810 


138 


-1.9811 


surface 


642 


-1.5617 


560 


-1.5425 


633 


-1.6251 


627 


-1.6181 


bulk 


1415 


-2.0120 


1415 


-2.0074 


1424 


-2.0040 


1446 


-2.0038 


total 


2057 


-1.8714 


1975 


-1.8755 


2057 


-1.8874 


2073 


-1.8871 



14 i- 
13 - 
12 - 
11 - 




l. ■ ■ 



-1.888 -1.887 -1.886 -1.885 -1.884 -1.883 -1.882 

Noncrystalline Crystalline 

<E> (eV/atom) 

FIG. 3: Histogram showing the distribution of energies ((E)) 
for 25 resolidified 2057 atom clusters. Two resolidified clus- 
ters were identified as fee truncated octahedra (with ener- 
gies of approximately -1.8820 eV/atom) while the remaining 
23 were identified as icosahedron-like (energies below -1.8840 
eV/atom). 

number of extra {100}-facets (see figure |]) distributed 
about the surface. Consequently, the binding energy of 
the surface atoms of the resolidified icosahedron is sub- 
stantially larger than that in both the Mackay and anti- 




FIG. 4: A 2057-atom resolidified icosahedra. Note the anti- 
Mackay type surface reconstruction but also the extra {100}- 
facet. This structure had the highest binding energy of any 
cluster produced in the resolidification trials. 



Mackay icosahedra. 

A closer examination of the resolidified icosahedra re- 
veals how they differ from anti-Mackay icosahedra. Re- 
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B 





FIG. 5: Tetrahedron with a regular Mackay surface termina- 
tion (A), an anti-Mackay surface termination with a stacking 
fault in the penultimate layer (B), and tetrahedron with a 
twin plane one layer lower (C). 



call that a Mackay icosahedron can be constructed from 
twenty fee tetrahedra. Likewise, an anti-Mackay icosa- 
hedron can be constructed from twenty tetrahedra, each 
of which has a surface with atoms in hep positions, as 
shown in figure pi We refer to the tetrahedra that make 
up a Mackay icosahedron as type A. The tetrahedra that 
make up an anti-Mackay icosahedron will be referred to 
as type B. Note that the type B tetrahedra has a stack- 
ing fault in the penultimate layer (...ABCABA) as the 
surface atoms lie in hep positions. The resolidified icosa- 
hedra have been found to consist of a mixture of type 
B tetrahedra, and a third type of tetrahedra, which we 
will refer to as type C. Type C tetrahedra have a twin 
plane in the third shell from the surface (...ABCACB). 
This third type of tetrahedra is also shown in figure ^|. 

The arrangement of these type B and type C tetrahe- 
dra for the 2057-atom resolidifed cluster from figure ||, 
is shown in figure ^|. The extra {100}-facets, visible in 
figure [|, occur at some of the edges between type B and 
type C tetrahedra. Table || compares the energies of the 
2057-atom resolidified icosahedron, cuboctahedron and 
Mackay icosahedron. It is clear that this surface recon- 
struction considerably lowers the surface energy of the 
resolidified icosahedra. Overall, the 2057-atom resolid- 
ified icosahedron has a total energy per atom that is 5 
meV lower than the cuboctahedron, and 16 meV lower 
than the Mackay icosahedron. 

Table III further decomposes the total energy of 2057- 
atom structures into the glue and pair potential compo- 
nents. The resolidified icosahedron is able to achieve a 
considerably lower glue energy than the two conventional 
structures, which more than compensates for an increase 
in the pair energy. We can further decompose the pair 
energy into two partsa: 



e + E stTa 



(5) 



where n nn is the number of nearest neighbors, e is the 
depth of the pair potential and -Entrain is the energetic 
penalty for pair distances that deviate from r m j n , the 
pos itio n of the minimum of the pair potential </>. Ta- 
ble III shows that the increase in the pair energy of the 
resolidified icosahedra comes chiefly from an increase in 
the strain energy. Thus, the surface reconstruction of the 
resolidified icosahedra is able to considerably improve the 




FIG. 6: This figure shows two arrangements of type-B and 
type-C tetrahedra considered here. The top pattern (with 
m = 7 type-C tetrahedra) shows the arrangement of tetra- 
hedra in the 2057-atom resolidified icosahedron shown in fig- 
ure ^. The lower pattern (with m = 12 type-C tetrahedra) 
shows the arrangement of tetrahedra in the icosahedra which 
appear to be stable at larger sizes. 



surface energy via the glue term, incurring a smaller in- 
crease in strain energy. 

We now wish to examine how the surface reconstruc- 
tion improves the glue energy. It is instructive to look 
at the culmulative contribution to (n^), from pairs with 
Tij < r for atoms i on the surface (defined here to be 



atoms with ni < 11): 



(nf(r)) ,= 

\ % K •'/surf 



surf 

£ 



p{r lj ). 



(6) 



This quantity, and the corresponding glue energy, are 
compared in figure [?| for the 2057-atom cuboctahedron, 
Mackay icosahedron and resolidified icosahedron (note 
that the energy curve follows the shape of the {nf (r)) sur{ - 
curve as U is approximately linear away from its min- 
imum). After the contribution to (rij) from the first 
shell (jij < 4.25 A say), the Mackay icosahedron has 
the largest (nf) f value. This is to be expected as the 
Mackay icosahedron has only {111} facets at the surface. 
However, for the resolidified icosahedron (n < (r)) c be- 

' \ % V JI surf 

comes largest beyond r = 4.6 A where the contribution of 
the second-nearest neighbors begins to make an impact. 
Note that since n^(r cu t) = Jij, it is clear from the figure 
that the resolidified icosahedra will have the lowest glue 
energy at the surface. 

The relatively large contribution from next-nearest 
neighbors to rii, arises due to the small differeiiee in en- 
ergy between the {111} and {100} faces for leadliJ. While 
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TABLE III: A comparison of the surface, bulk and total energies per atom of a 2057-atom cuboctahedron, a 2057-atom closed- 
shell anti-Mackay icosahedron and the 2057-atom resolidified icosahedron in figure Ell Also shown is the decomposition of the 
total energy into -E" pa ir + -Egiue (equation and the strain energy S s tr a m (equationjEJ) . 



Cluster 


(-^bulk) 


(-^surface) 


(E) 


(%ue) 




{-^strain) 










(eV/atom) 








Cuboctahedron 


-2.0208 


-1.5901 


-1.8817 


-1.7368 


-0.14492 


0.01609 


11040 


Mackay Ico 


-2.0141 


-1.5510 


-1.8714 


-1.7253 


-0.14614 


0.01802 


11256 


Resolidified Ico 


-2.0032 


-1.6275 


-1.8876 


-1.7522 


-0.13544 


0.02634 


11093 



Resolidifed 




r(A) 



FIG. 7: A comparison of (nf (r)) f and U ((nf(r)) f ) for 
the m = 7 2057-atom resolidifed icosahedron (solidj in fig- 
ure ^, the 2057-atom cuboctahedron (dashed) and the 2057- 
atom Mackay icosahedron (dash-dot). The cut-off, at which 
nf{r C Mt) = n i, is indicated by the vertical dashed line. 



an atom in a {100} face has fewer nearest- neighbors, 
it has more next-nearest neighbors than an atom in a 
{111} face. The glue energy reflects this small differ- 
ence in energy by weighting the contributions from next- 
nearest neighbors more heavily than might otherwise be 
the caseo. This would seem to be the key feature of the 
potential that stabilizes the novel surface reconstructions 
of the resolidified icosahedra. 

In the next section, we will explicitly construct new 
icosahedra from type B and C tetrahedra. This will en- 
able us to compare these new icosahedra to fee structures 
at cluster sizes where it is currently too expensive to con- 
duct repeated resolidification trials. 



V. NEW ICOSAHEDRA 

Using the two types of tetrahedra identified in the 
resolidified icosahedra, we can construct 2 20 icosahedra. 



However, many of the 2 20 possible icosahedra can be 
identified after a rotation. In the appendix , we show 
that there are only 17284 unique ways of constructing an 
icosahedron from the two types of tetrahedron. 

This a large configuration space to search for the best 
arrangement of tetrahedra. However, a pair of clusters 
which arc mirror images of one another will be energet- 
ically equivalent i.e. the pair will be chiral isomers. In 
fact, there are 1048 of the 17284 clusters which are invari- 
ant under reflections, leaving 8488 pairs of chiral isomers. 
This reduces the number of energetically distinct clusters 
to at most 9536. Furthetr-jthe resolidification trials de- 
tailed in Hcndy and Hallia are a way of sampling this 
configuration space to discover low-energy clusters. 

To construct one of these new icosahedra from scratch, 
we begin with an (n — 2)-shell Mackay icosahedra. At 
faces where type-B tetrahedra are desired, add a fur- 
ther (n — l)-shell Mackay icosahedral face, and at faces 
where type-C tetrahedra are desired, add an (n — 1)- 
shell anti-Mackay face. To complete the penultimate 
shell, Mackay icosahedral edges and vertices are added. 
Now to complete the outer shell, add an n-shell anti- 
Mackay face at faces where type-B tetrahedra are desired 
(the corresponding type-B tetrahedra are now layered 
as . . . ABCABA), and an (n + l)-shell anti-Mackay face 
where type-C tetrahedra are desired (the corresponding 
type-C tetrahedra are now layered as . . . ABCACB). 

The number of atoms in a n-shell icosahedra (n > 2) 
constructed as above from m type-C tetrahedra and 20 — 
m type-B tetrahedra is then given by: 

NewIco(n, m) = Ico(n — 2) 

+ 30n - 48 + m(n 2 - n + 1) 

+ (20-m)(n 2 -3n + 3) (7) 

The cluster from figure ^, which was the best structure 
produced in the trials shown in figure |3|, has seven type-C 
tetrahedra and thirteen type-B tetrahedra. With m = 7, 
equation ^ gives a sequence of 309, 565, 931, 1427, 2073, 
2889, 3895, 5111, 8253 ... atoms in each closed-shell 
cluster. We note that a n-layer type-C tetrahedra con- 
tains 2(n — 1) more atoms than the type-B tetrahedra 
and, hence a n-shell icosahedron with m type-C tetrahe- 
dra has 2m(n — 1) more atoms on the surface than an n- 
shell anti-Mackay icosahedron (which is constructed from 
20 type-B tetrahedra). Thus, by increasing the propor- 
tion of type-C tetrahedra, the density of surface atoms is 
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FIG. 8: Energies of clusters versus size relative to the fit to 
the cuboctahedra sequence: m — 7 icosahedra (open gradi- 
ent symbol), m = 12 icosahedra (open box) and truncated 
octahedra (filled circle). 



increased, relative to the anti-Mackay icosahedra, as seen 
in table M. 

In table 0, the energetics of a 2073-atom to = 7 clus- 
ter, constructed using the design in figure ^ and then 
relaxed, can be compared to those of a 2057-atom reso- 
lidified icosahedron. From the CNA-classification of the 
atoms in these clusters, it is clear that the two clusters 
are very similar in structure and energetics. Although, 
we did not explicitly include the extra {100}-atoms at 
the outer edges in the above construction, which appear 
as CNA signature C in the table, these arise as the con- 
structed cluster relaxes. An interesting feature of the 
structure of this cluster is that it has a chiral isomer, as 
the m = 7 pattern ^ does not have a mirror symmetry. 

Figure ^ shows the energies for a sequence of m = 7 
icosahedra, relative to the cuboctahedra sequence, with 
the first design given in figure ^. From figure ^| it appears 
that this type of icosahedra is the lowest energy structure 
for sizes from 900-5000. The largest cluster of this design 
that lies above the interpolated fit to the truncated oc- 
tahedra sequence is the 5111-atom cluster. 

We have found one other design which appears to be 
stable at larger sizes. This design is also shown in figure ^ 
and has 12 type-C tetrahedra (to = 12). The energies of 
this sequence are shown relative to the cuboctahedron fit 
in figure @. The first energetically favored m = 12 icosa- 
hedron occurs at a size of 5211 atoms. This sequence 
continues to give the most energetically favored clusters 
we have found until a cluster size of 18097 atoms, where 
the truncated octahedron sequence appears to become 
more stable. Thus, this sequence gives the most ener- 
getically favored structures known from a size of 5211 to 



15191 atoms. 



VI. DISCUSSION 

We have described a new family of icosahedral struc- 
tures, which include the anti-Mackay icosahedra. The 
family is characterized by a Mackay icosahedral core and 
a surface reconstruction that can extend into the two 
outer layers. They can be explicitly constructed from two 
types of fee tetrahedra: the first type of tetrahedra has an 
anti-Mackay surface termination (type-B), and the sec- 
ond type of tetrahedra has a twin-plane fault in the third- 
outermost layer (type-C). The anti-Mackay icosahedra is 
constructed entirely from type-B tetrahedra. Structures 
in this family of icosahedra were shown to be the lowest 
energy structures known for the lead glue potential ([!]) for 
certain sizes. This study illustrates how new structural 
forms can arise when many-body effects are included in 
the interatomic potential. 

Icosahedra constructed with a sufficient number of 
type-C tetrahedra possess a density of surface atoms 
comparable to that of the Mackay icosahedra, but with- 
out the low coordination atoms at the edges and vertices 
of a Mackay icosahedron. Such a surface reconstruction 
was shown to be favored by the glue-term in the inter- 
atomic potential (jl]) over fee structures largely due to the 
effects of next-nearest neighbors. We commented that 
the relatively large effect that the next-nearest neighbors 
have on the glue term is due to the small difference in 
energies between the {111} and {100} faces. This next- 
nearest neighbor contribution is also the feature of the 
potential that favors fee structures over Mackay icosahe- 
dra at all sizes, as noted by Lim, Ong and ErcolesslL-a in 
their original study of lead clusters. 

While we have not exhaustively searched all possible 
icosahedral structures that are part of this new family, we 
have identified two configurations of type-B and type-C 
tetrahedra that show particularly favorable energetics. 
The first structure (which occurs in a pair of chiral iso- 
mers) has seven type-C tetrahedra and is the lowest en- 
ergy structure known for the potential (Q) over a size 
range of 900-5000 atoms. The second structure (which is 
symmetric under reflection) has twelve type-C tetrahedra 
and is the lowest energy structure known for the poten- 
tial (|l|) over a size range of 5000-18000 atoms. Thus, 
above 900 atoms, fee structures are not favoured by (|l|) 
until at least cluster sizes of 18097 atoms, where the trun- 
cated octahedra appears to be favored. This is probably 
a conservative lower bound on the size where large fee 
structures appear as global minima of the potential ([!]) . 

These results continue to emphasize the tendency of 
this potential to produce non-fee structural forms. A 
recent global minimization of this potential for cluster 
sizes of up to 160 atoms, found that these clusters do not 
adopt fee structures at any size in this rangeEM. Further, 
it seems likely that no globally minimum fee structure ap- 
pears in the window between that study and the results 
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presented here (i.e. between 160 and 900 atoms) although 
we have not explicitly jiemonstrated this here. Simula- 
tions of lead nanowiresa have also shown the emergence 
of non-fcc structures. In all these cases, the non-fcc char- 
acter is related to the small difference in energy between 
the {111} and {100} faces (although, disturbingly, it ap- 
pears that the cut-off distance r for the. glue energy also 
plays a role, at least for small clusterstiJ). 

We have focused here on icosahedral structures because 
they were produced in the melting and freezing simula- 
tions of Hendy and HallllS. While these structures were 
found to be lower in energy than other known structures, 
there is no guarantee that they are globally minimum. 
Indeed, simulations of freezing often produce icosahe- 
dral structures preferentially, either for thermodynamic 
or kinetic reasons, irrespective of whether these struc- 
tures ate globally optimal (see for example Chushak and 
BartellEa). It is possible that the novel surface recon- 
structions seen here may also stabilize decahedral forms, 
for example, which have not been seen in freezing sim- 
ulations, but which may compete energetically with the 
icosahedral structures seen here. 



APPENDIX: N-COLORED ICOSAHEDRA 



Thus, we can determine C by determining \X\ g for each 
g £ G. In fact, for each jeG, \X\ g = N m , where to < 20 
is the number of orbits of faces under g. Determining the 
value of to for each g is straightforward, and can be done 
using a cardboard cut-out icosahedron (or counted on a 
computer). 

Doing so, we arrive at the formula for C for the N- 
colored icosahedron: 

C = — (TV 20 + 247V 4 + 207V 8 + 15TV 10 ) (A.2) 
60 

Thus for TV — 2, the number of 2-colored icosahedra is 
17824. 

We can also consider N-colored icosahedra under re- 
flections as well as rotations. The size of the symmetry 
group is now doubled as we add a generator of reflections 
to the group: \G\ = 120. The formula for the number of 
N-colored icosahedra unique under rotations and reflec- 
tions is: 

C = (jv 20 + N w + 24(^2 + 7v 4 ) 

+ 20(TV 8 + TV 4 ) + 15(TV 12 + TV 10 )) (A.3) 



We consider painting an icosahedron using up to N 
colors, where each face can only be painted a single color. 
There are clearly TV 20 ways to paint the icosahedron, but 
some of these painted icosahedra will simply be rotations 
of other painted icosahedra. 

To count the numberp-af unique painted icosahedra we 
use Burnside's theoremEa. Let X be the set of colorings 
of an icosahedron (\X\ = TV 20 ) and G be the rotational 
symmetry group of the icosahedron ( \G\ = 60, consisting 
of the identity, 15 180° rotations about edges, 20 120° ro- 
tations about faces, and 24 72° rotations about vertices) . 
Now for each g 6 G, we define X g — {x e X \ gx = x}, 
where gx is the new coloring obtained by rotating the 
coloring x via the rotation g. 

Now, Burnside's theorem gives the number of color- 
ings, C, that are unique when acted upon by the finite 
group G: 



(A.l) 



Thus, if we identify mirror images, formula (A.3) applies, 
and hence the number of 2-colored icosahedra is 9436 
(TV = 2). Of the 17824 2-colored icosahedra, 1048 are 
invariant under reflection symmetry, while the remaining 
16776 2-colored icosahedra come in 8388 pairs of mirror 
images. 
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